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Two Dimensions. 
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We study the Principal Chiral Ginzburg-Landau- Wilson model around two dimensions within the 
Local Potential Approximation of an Exact Renormalization Group equation. This model, relevant 
for the long distance physics of classical frustrated spin systems, exhibits a fixed point of the same 
universality class that the corresponding Non-Linear Sigma model. This allows to shed light on the 
long-standing discrepancy between the different perturbative approaches of frustrated spin systems. 
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There is now a general agreement about the field theo- 
retical treatment of the SO(N) spin system. The pertur- 
bative approaches performed around four dimensions on 
the Ginzburg-Landau- Wilson (GLW) model, around two 
dimensions on the Non-Linear Sigma (NLcr) model and 
in a 1 /N expansion give a consistent picture of the criti- 
cal physics of this system everywhere between D = 2 and 
D = 4 H. This picture has also been confirmed by non- 
perturbative methods based on truncations oLWilson Ex- 
act Renormalization Group (RG) equations 0[J. Amaz- 
ingly, there is no such agreement for many systems whose 
symmetry breaking pattern is not given by S.QIN) — > 
SO(N — 1) among which.-are-.superfluid 3 He Erfc 
trated antiferromagnets O cjL.— superconductors li^ 
electroweak phase transition etc. Generically per- 

turbation theories predict that these systems undergo a 
first order phase transition npar D j= 4 and a second or- 
der one around D = 2 HByMtJ. The origin of this 
discrepancy is not yet understood and calls for a non- 
pcrturbative approach. 

In this paper we study, by means of the Wilson Renor- 
malization Group approach, the Principal Chiral (PC) 
model, corresponding to the symmetry breaking scheme 
SO{N) <g> SO(N) -> SO(N), which is the simplest one 
exhibiting the non trivial features previously quoted. 
The PC model is the low energy effective field theory 
of a whole class of systems among which frustrated an- 
tiferromagnets. A particularly important example is 
the Heisenberg antiferromagnet on the triangular lattice 
(AFT). Due to the frustration, the order parameter is a 
triad of orthonormaJ-Ypctors, i.e. a 5*0(3) rotation ma- 
trix R = (ci, 62, 63) E3'E-3 . We consider, in the following, 
the generalization to N orthonormal vectors e^'s with 
N components, i.e. SO(N) matrices. The long distance 
physics of this generalized AFT model is thus equivalent 
to that of orthonormal frames interacting ferromagneti- 
cally: 
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<i,j> a=l 



<i,3> 



The Hamiltonian (|]J) is invariant under the SO(N) (g) 
SO(N) group of left U 6 SO(N) and right V £ SO(N) 
global transformations: R l — > UR l V. Since, in the low 
temperature phase, the residual symmetry group con- 
sists in a (diagonal) SO(N), Eq. (|l|) is indeed a lat- 
tice version of the PC model. Whereas the microscopic 
derivation for frustrated spin systems leads in general to 
anisotropic interactions between the e a 's, i.e. J is a- 
dependent, we consider here the isotropic case where all 
the J a 's are equal. It has been shown for a large class 
of frustrated spin systems, among which the AFT model, 
that the anisotropy is anyway irrelevant, at least near two 
dimensions, for the critical properties we are interested 
in IM 

Let us first sketch out the experimental and numerical 
situation for frustrated spin systems which, in D = 3, is 
far from being clear. Indeed, the behaviour of systems 
that are supposed to be described by the PC model like 
AFT (CsVCl 3 , RbNiCl 3 ) and Helimagnets (Ho, Dy, Tb) 
are affected by the presence of disorder localized near 
the sample surface and, possibly, by topological defects. 
As a consequence, the critical eptppiients strongly vary 
from one compound to another OEJ. Numerically, the 
situation is also confused since simulations performed on 
the PC model and directly on the AFT model lead re- 
spectively to first order and second order behaviour with 
exponents of an unknown universality class E3 . 

Beyond this apparent lack of universality at the ex- 
perimental and numerical level, the theoretical situation 
already exhibits the puzzling features previously men- 
tioned. Around D = 2, the critical physics is obtained by 
means of a low temperature expansion performed on the 
PC NLcr model. A fixed point is found for any N > 2 in 
D = 2 + e dimensions so that a second order phase tran- 
sition is expected ffl. On the other hand, the weak 
coupling expansion performed in D = 4 — e on the PC 
GLW model suggests a first order phase transition since 
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no fixed point is found for any N > 2 U. As such, 
the situation is not paradoxical since perturbation theo- 
ries are only trustable in the immediate vicinity of their 
respective critical dimension. However if, as usual, we 
extrapolate the perturbative results to D = 3, the two 
results come into conflict. It is thus of first importance to 
clarify this theoretical situation before hoping to describe 
real materials. 

From the theoretical point of view, the crucial fact is 
that the calculation of the f3 functions in the two different 
perturbative approaches relies on qualitatively different 
grounds. Indeed, the (3 function of a NLa model built on 
a manifold G/H i©nly depends on the symmetry breaking 
scheme G — > H & - i.e. on Goldstone modes - whereas 
that of the GLW near D = 4 is sensitive to the repre- 
sentation of G spanned by the order parameter chosen to 
realize the symmetry breaking scheme. This feature can 
be fully appreciated in the N — 3 PC model. Indeed, 
since 5*0(3) ® 50(3) is isomorphic to 5*0(4) the symme- 
try breaking pattern is that of the usual four component 
spin system: 50(4) — > 50(3). The perturbative f3 func- 
tion of the A = 3 PC NLa model in D = 2 + e is thus 
identical to that of the A — 4 vector model, although 
the symmetry breaking scheme is realized with a rotation 
matrix which is_a-SO(4) tensor and not with a four com- 
ponent vector EMLJ . If this perturbative result remained 
true beyond D = 2 + e, as it is believed in the SO(N) 
vectorial case, we could expect the critical behaviour of 
the PC model to be determined by the same fixed point 
as the A — 4 vector model everywhere between two and 
four dimensions. This is however not the case, at least 
perturbatively in the vicinity of D — 4, since there is no 
fixed point in the GLW approach. 

The origin of the discrepancy between the two ap- 
proaches can be ultimately traced back to the (non per- 
turbative) spectrum of both models. Whereas it is very 
likely that in the SO(N) case with a vectorial order pa- 
rameter the NLcr and GLW models share the same low 
energy degrees of freedom everywhere between two and 
four dimensions, it is no longer the case for models with 
more general order parameters and symmetries. For ex- 
ample, for the N = 3 PC model, the spectrum of thf 
D = 2 NLcr model consists in four massive particles E3 
whereas the spectrum of the D = 4 GLW in the high 
temperature phase involves nine massive particles. Ide- 
ally, we should understand at a non-perturbative level 
how these two field contents are linked together in D = 3 
and how they are related to the degrees of freedom of the 
underlying microscopical system. This is a formidable 
task that will not be tackled here. 

The question we address here is the possibility of a 
matching between the NLcr and GLW approaches when 
varying the dimension. This allows, at the same time, to 
test the validity of the NLcr model for frustrated systems, 
at least around D = 2. Indeed, due to the discrepancy 
between the two perturbative approaches and the ab- 
sence of experimental and numerical evidence of an 0(4) 
critical behaviour, the reliability of the NLa model ap- 



proach has been questionncd @ . Clearly, the answer to 
these questions escapes a perturbative treatment. In gen- 
eral, the 1/N expansion provides a powerful tool to link 
up different perturbative methods. In the case of matrix 
models such an analysis is however plagued by technical 
difficulties. Some progress have been recently obtained 
but arc confined to the leading order EaH. The Wil- 
son's scheme, .wjhich has been successfully used in many 
topics E^I^EjO'EI, turns out to be the most efficient ap- 
proach. In this paper, we study the PC GLW model near 
D = 2 by means of the Wilson - Polchinski Exact Renor- 
malization Group within the Local Potential Approxima- 
tion (LPA). We show that the GLW and NLcr approaches 
can be reconciled in the vicinity of two dimensions. More 
precisely we show by a RG analysis that the two models 
belong to the same universality class near two dimensions 
since the GLW model exhibits a non trivial fixed point 
identical to that of the NLa model. 

The partition function of the PC GLW model is ob- 
tained by writing the most general SO(N) (g) SO(N) in- 
variant potential, at most quartic in A x A real matrices 
M, that favours orthogonal matrices for the field M: 



Z = / DM exp 



J d D x ^Tr (V*M.VM) 



+ -Tr 'MM + fiTr ('MM) 2 + A(Tr 'MM) 2 



(2) 



The domain of parameters of interest for us is given by 
A > since, in this case, the minimum of the poten- 
tial in the broken phase is given by M(x) = Rq where 
i?o £ 50(A). In this phase, the model displays a 50(A) 
symmetry, so that the symmetry breaking scheme is 
50(A) ® 50(A) -> 50(A) and thus corresponds to the 
GLW version of the PC model. 

Our aim being to make contact with the NLa model, 
let us show how the orthogonality of the lattice order pa- 
rameter of ([j]) can be recovered from (||). Let r and \x 
go to infinity, the ratio r/4fi being fixed. In this limit, 
one recovers the partition function of the PC NLa model 
where a delta function enforces the orthogonality con- 
straint on M at each point: 



Z = J DM cxp-i f d l 



Tr (V'M.VM) 



(3) 



-2/xTr ( 'MM + — ) + 2A(Tr 'MM) 2 



-> J DM S (*MM - 4^ ex P~^ J d ° x Tr (V'M.VM) 

(4) 

up to an overall constant. The quantity l/.g 2 = — r/4/x 
which corresponds to the minimum of (^) (when A <C /i) 
identifies with the inverse temperature of the NLa model. 
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Of course, since the preceding limit is made on the bare 
action, it does not allow to conclude how both models 
are related under RG transformations. We shall show 
that, around two dimensions, the GLW and NLer models 
actually converge to the same renormalized trajectory in 
the continuum limit. 

To realize this program we now study the evolution of 
the PC GLW model under RG transformations within 
the LPA. This approximation consists in truncating the 
effective Wilsonian action to its potential part V(M) = 
J d D x v(M(x)). Note that the LPA thus misses the field 
renormalization. The Wilson-PolchiHski equation for the 
potential density v(M) is given by Ha: 



dv „ 
di = Dv 



D-2 , 1 „ 

2 M ij V ij + ~^ V ij,ij 



(5) 



where = 5v/5Mij and t = In A, A being the dimen- 
sionless running scale. 

There are two different ways to exploit Eq. (H|). The 
first one is to search for an exact solution in any dimen- 
sion, having recourse to numerical integration. This pro- 
vides a powerful way to obtain precise values for critical 
quantities 0. The second one is to solve Eq. (|^) in a 
low temperature expansion. We follow this route since 
we are interested in qualitative features of the RG flow 
and we want to make contact with the standard pertur- 
bative analysis of the NLcr model. Mitter and Ramadas 
used the same techniques in the SO(N) case for a Pfopf 
of perturbative renormalizability of the NLcr model E3 . 

Let us parametrize the potential density v by v(M) = 
u (x)/9t w hh x = 9t t MM. In a perturbative approach 
the running potential has always a minimum as it is the 
case for the initial potential in (|) for t MM = l/g%. The 
running temperature <? t is thus defined via: 



du 
dx 



0. 



X = l 



(6) 



We now write the Wilson-Polchinski equation for the po- 
tential density u within the LPA: 

Oil 

— =Du-(D-2) Xij u' ji 

~ ( u 'jl u 'jk + u 'ji u 'kj + u 'ij u jk + u 'lj u 'kj) Xlk 

[( u %.jp + u "j,jp + u js, P j + u "j,pj) Xs P + 2iVi4 



+9t 



dt 



9t 



(7) 



Under the iterations of the RG, all SO(N) <g> SO(N) in- 
variant terms are generated so that the evolved potential 
writes: 



U (X {^pi,9i;. ••;p«,g n (01) — \>i,«i;...;Pi,«(*) 

» {Pk,Qk} 

[Tr ( X -l)^r...[Tr . 

(8) 

The Wilson-Polchinski equation (^) generates the flow of 
all the ^p lt qi;...;p„,q n (tys. When combined with (JsJ) we 
also get the evolution of g t : 



d -§ = -{D-2)gl 



9t 



4tt 2A 2 ,i(t) + 2AAi j2 (i) 



(12N + 12)A 3 ,i(t) + 24AfA lj3 (t) + 4(A 2 + N + 4)A a ,i;i,i(t) 



+4(2N + l)A 2 ,i(<) + 4(A^ 2 + 2)Ai >2 (t) 



(9) 



The flow analysis shows that all the X pi , qv .,. ;Pn 9n (i)'s 
are irrelevant coupling constants: after an exponentially 
rapid transient regime, their scale dependence is entirely 
controlled by that of gt: 

^Pi,qi;—;p n ,q n (t) ~* ^pi,?i;-» ;p»,«n ^pi,8i;— vp n ,in 9t + ^(fft) 

(10) 

Therefore, for any initial conditions, the flow is 
driven towards a one-dimensional renormalized trajec- 
tory parametrized by gt whose evolution, obtained from 
(ft) and (HO), is given at leading orders by: 



N-l 

47T 



9t + 0(gt) 



(11) 



This P function identifies with that of the temperature 
of the PC NLcr model calculated perturbatively O . It 
however differs from the standard expression where the 
coefficient N — 1 is replaced by N — 2. The origin of this 
difference is that, within the LPA, the field renormaliza- 
tion is set equal to one. If the field renormalization had 
been taken into account, which is the case in the next or- 
ders in the derivative expansion, we would have obtained 
the correct coefficient. This difference is irrelevant for 
our purpose. 

Let us indicate how, in two dimensions, our previous 
results allow to recover, in the continuum limit, the hard 
constraint of the NLcr model (f|). After the transient 
regime - i.e. Eq. ([[(]) - u(x, {^pi,9i;...;p„,«„(*)}) nas con- 
verged towards u(x, gt) which can be expanded in powers 
ofg?: 



u(x,9t) = Y,(9t) k n {k Hx) 



(12) 



k>0 



We have found the exact form of (x) so that the dom- 
inant part of the potential density at low temperature 
writes: 
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v(M) 



07 



? 'MM - 1 



(13) 



Suppose now that, after blocking, the model having con- 
verged to the one-dimensional renormalized trajectory, 
the effective running temperature has reached the value 
g^, at scale fi. Reversing the flow on this trajectory, 
towards the ultraviolet, Eq. (|Tl]) gives the bare temper- 
ature g at scale of the overall cut-off A (typically, the 
inverse lattice spacing). Due to asymptotic freedom go 
goes to zero when taking the continuum limit Ao — > oo. 
It is easy to see from ( |l3| ) that, in this limit, the configu- 
rations contributing to the partition function correspond 
to SO(N) matrices (up to a normalization): the delta 
constraint of (Q) is recovered from RG transformations. 
Thus, in the continuum limit, the GLW and NLcr models 
coincide. The statistical interpretation of this is that the 
soft field GLW model appears as the block-spin iterated 
NLcr model. 

These results show that the PC GLW and NLcr mod- 
els belong to the same universality class near two dimen- 
sions. This is a strong evidence of the validity of the NLcr 
model approach and of the existence of a second order 
phase transition near two dimensions. Thus, the criti- 
cal behaviour of the PC GLW model must change as D 
varies between D = 2 and D = 4. This, of course, relies 
on the assumptions that our results persist beyond the 
low-temperature expansion and the LPA, and that the 
e = 4 — D expansion of the GLW model is meaningful. 
The change of critical behaviour could be a general fea- 
ture of models that are afflicted by the same troubles even 
if their origins - presence of topological excitations, role 
of irrelevant operators - certainly depend on the precise 
model under study. In any case, analyzing this requires to 
vary the dimension and to use the nest ocders of approx- 
imation in the derivative expansion H~H . A somewhat 
similar study has been performed for superconductpis 123 
and for the Kosterlitz-Thouless phase transition CJ . In 
the context of the PC model, it will be addressed in a 
future publication. 
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manuscript. 
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